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A full order adaptive observer is described 
for observing the states of a single-input single- 
output observable continuous differential system 
with unknown parameters. Convergence of the 
observer states to those of the system is accom- 
plished by directly changing the parameters of 
the observer using an adaptive law based upon 
Liapunov stability theory. Observer eigenvalues 
may be freely chosen. Some restriction is placed 
upon the system input in that it must be suffi- 
ciently rich in frequencies in order to Insure 
convergence. 


Introduction 

The Luenberger observer [1,2,3] allows 
extraction of all the states of an observable 
linear system when given the output and the 
parameters of the system. In some cases the 
system parameters may not be known; consequently, 
in these cases the state observations are sub- 
ject to error. Previous investigators of this 
phenomenon have attempted to estimate the error 
[4] or change the observer parameters in some 
beneficial way [5]. Their analysis suffers in 
that the error cannot be guaranteed to vanish. 

We report a full order observer for a restrictive 
class of systems (that is, single-input single- 
output observable continuous linear differential 
systems in the absence of a deterministic or 
random disturbance vector) for which the observa- 
tion error is guaranteed to vanish regardless of 
the size of the constant or slowly varying para- 
meter ignorance. The observer parameters are 
directly changed in a way that satisfies a quad- 
ratic Liapunov function of the error and the 
correct but unknown Luenberger observer para- 
meters. The observer poles may be placed at any 
stable location and no derivatives are required 
in the adaptive law. 


some or all of the elements of matrices A and B 
are unknown, _A is stable, w® may be unknown, and 
the pair (C,A) is completely observable. The 
observer is of the form 

z = Fz + GCw + Dr + Hu z(0) = z^ 

( 2 ) 

F nxn G nxl 

D nxl H nxn and diagonal 

where z is arbitrary and u is a control vector 
yet to be defined but with the property that u-*-0 
as t-«». The problem is to adaptively form a 
triple (G,DjT) so that the error vector defined 
as e = z-T w vanishes as the system adapts. T 
is a non-singular square matrix with the property 
that CT=C. 

Define a transformation x •= T ^w so that e =z-x. 
Then (1) becomes 

x = A Q x + T _1 Br x(0) = Tw° 

y = CTx “ Cx (1A) 

A Q - T _1 AT 
and (2) becomes 


z = Fz + GCx + Dr + Hu 
z(0) « Z° 


(2A) 


It is desired that A. •= T *AT be in the "output" 
form 



The Problem 


A differential system is assumed of the 


form 


w = Aw + Br 
y - [10 0 0]w 


w(0) = w 


( 1 ) 


A nxn 
B nxl 


for which only the single output y= Cw = w^ is 
available for measurement . It is assumed that 


-a nl 0 0 0 ... Oj 

wherein the first column contains the system para- 
meters and all other elements are zero save the 
super diagonal elements, which are unity. It is 
clear that for any non-zero matrix A there corres- 
ponds a similar matrix A^, although the elements 
of the similarity transformation may be unknown 
if elements of A are unknown. The following 
theorem defines the additional restriction that 
must be placed upon A so that both Aq = T - *AT and 



Aq-GC, a constant matrix 


CT = C. 

Theorem [proof given in ref. 11J . Let A be 
an nxn matrix, C = [k,0,0, . . . ,0] a lxn matrix 
with k^O, A^ an nxn matrix in output form, and 

T = {. . . ] an nxn nonsingular matrix. There 
T ■ 

exists an (n-l)xn matrix T such that A=TA^T iff 
the pair (C,A) is completely observable. As a 
result of the theorem, any observable system (1) 
may be placed by similarity transformation into 
system (1A) with CT=C. The elements of T may be 
unknown since A is unknown. The problem will 
be considered as defined by equation (1A) and 
(2A), so that e = z-x must vanish. Eventually 
the problem of constructing w from x will be 
solved . 


The Adaptive Law 

It is now assumed without restriction that 
some stable "nominal" plant matrix is known or 
is chosen so that A^ ° A^ + AA^ where has all 

known elements and is in output form. Then AAq 

contains all zero elements except for the first 
column which has unknown elements. The vector 
error equation may then be written, letting 
e = z-x, as 

e » Fe + (F + GC - A Q - AA Q )x + ABr + Hu 

where AB “ D-T X B. F may be chosen as F “A^-GC . 
The resulting error equation is 

e = (A Q -GC)e - AA Q x + ABr + Hu (3) 

It is desired to reduce AA^, AB, and u to zero. 
Then, if G is chosen so that A^-CC Is stable, 

e(t) approaches zero. A theorem of Luenberger 
[1] allows the eigenvalues of A^-GC to be placed 
arbitrarily by selection of G, with the sole 
exception that A^-GC cannot have the same eigen- 
values as Aq. Consequently, A^-CC can always be 

made stable. Moreover, it is assumed throughout 
that G will be chosen so that all elements of 
Aq-GC will be constant even under changes in 

AAq due to the adaptive law. 

The error between plant states x and observer 
states z may be measured only by the scalar 
e^ = Zj-y = z i~ x i‘ T° insure that only avail- 
able measurements are included in the adaptive 
law, the vector error equation (3) is "collapsed" 
to yield a scalar differential equation of the 
form 


n-1 


V u U) V (i) 

l = L “i*, 

1=0 A i=0 


(4) 


i=0 


m . . n-1 

l V + l Vi 

i-0 


(i) 


where 

k i e K 0 - 

e -AAq and its several derivatives 

8^ e AB and its several derivatives 
h t t H 

Letting p = d/dt, the left side of (4) may be 
written as 


iSi (p + V e i 

Now a reduction of order technique, similar to 
that of Gilbart and Monopoli [6], will be applied. 
n-1 of the (p + A ) terms will be selected and 
factored out of the right side of (4) excluding 
the u^ terms. Assuming that p + A^ was not 

selected, the error equation then has the form 
n n n+n+1 

(p + A r ) ± n 2 (p+A^)e^ -^2 <P+A t ) [ I VJ 


-f <♦<» 


where 


n (p + a ) v 

j=2 3 1 


n-1 

)+ l 

i=0 


_<i> 


h i u i 


(i) 


(5) 


1=1,2, . . . ,n-l 


n (p + A )v 
j=2 3 


.(i-n) 


i«=n,n+l, . . . ,n+m 


^ are functions of {a^} and (B^l, and 
(k) 

v^ ) is a function of derivatives of but 
does not contain for any j . 


(i+1) 


CD 


J i ’ 


specifically, is made equal to the negative 

of all the terms in which appears for 

every i and j. Here it is noted that by construc- 
tion neither nor u^ require a derivative net- 
work for implementation. Then (5) becomes 


(p+V ± n 2 (p+A 1 )e 1 = ± n 2 (p+A^ 


• n-hn+1 

l 

i=0 


^i v i 


and each u. 


.(♦f\ v«) 


( 6 ) 


Taking Laplace transform of (6) and dividing by 


^(s + A^ yields 



Reconstruction of T 


t n+m+1 ■> 

Jo 


(7) 


+<initial conditlons)/^^ (s + X^) 
for which follows 

n+m-t-1 n 

e +1 e = l $ .v + I <[>. exp [-X.t], 
■ L i i=0 i=0 1 1 


Xj real 


Using the "nominal" matrix as initial condi- 
tions, the actual value of the system parameters 
may be determined by integrating the change in 
parameters {(Jk) until adaptation is complete and 
combining appropriately. This procedure may be 
accomplished while the adaptation progresses by 
forming a matrix l*(t) with elements composed of the 
combination of nominal values and the integrals 

/ C $.dt where {?.} is defined in (10). 
oi i 

u, the estimate of o>, is constructed from the 
observer output z by forming 3 1 z. Since lim z = x 

' V , - t -+<“ 

and lim T = T, so lim ui = to. 
t-K= 


where are constants depending upon the initial 

conditions if { X are distinct; otherwise some 

may be time dependent. (Note: should it be 

desired, when n is even, to have no real 
observer pole, the operation in (7) may still be 
made by modifying the right side of (7) in an 
obvious way) . 

^ 2 n+m+1 2 

A Liapunov function V = + ilj “i^l^ 

is chosen, and V is calculated. Following 
Shackcloth (7), V can be made to be of the form 

2 " 

v = - n^X.^ + e x [ exp l -X.t] (9) 

i“l 

ra o 

when <f> . v e for all i (10) 

i i 1 

Other adaptive laws can easily be chosen instead 
if it is desired to increase convergence speed 
[8.9]. 

From the form of V, e^ is stable in the 
sense of Lagrange with the region of attraction 
determined by the unknown constants and the 

exponential time function. Clearly the region 
of attraction shrinks exponentially with time 
and eventually vanishes; consequently e^ is 
eventually asymptotically stable, and lim e =0. 

£-+00 

All derivatives of e^ must vanish in the limit 
also since the error equation is linear of 
first order. Although the Liapunov function is 
defined on a non-compact manifold (i.e., con- 
tains e^ but not ft can be shown that 

{<>^} is eventually asymptotically stable if the 
input to the plant, r, is periodic and contains 
(n+m+l)/2 distinct frequencies, none of which 
has a phase shift of kit through the plant, k any 
integer. [10] It must be assumed that the adap- 
tive observer is limited to systems for which 
convergence to parameter differences is 

assured. 


Since £A.,£Bq,u^ approach zero, the vector 
e is eventually asymptotically stable if G is 
chosen so that A^-CC is asymptotically stable. 


Example 


A third order plant with one zero is con- 
sidered for illustration of the previous dis- 


] 

0 

i 

0 


r o 

V = 

0 

; ( W 

0 

-(•!**!) 

1 

-(a 2 +a 2 ) 

w+ 

c i 

C 2 


( 1 *) 


where C^, and are unknown. A 

transformation T that delivers the system into 
output form is 


.-1 


a 2 +a 2 


a l' ta l 


0 

1 


a 2+ a 2 


0 

0 

1 


Note that C = CT. Then in output form, (1*) 
becomes 



-(a 2 +a 2 ) 

i 

0 1 


0 

x « 


0 

1 

x + 

V S 1 


; (a o-V 

0 

0 


yv 


(1A*) 

y = X 1 = W 1 


The error equation is 



-(a 2 +g 2 ) 

1 

o' 


“2 

0 

o' 


’ 0 " 

e = 

-( a l+Rl) 

0 

1 

e + 

°1 

0 

0 

X + 

- 8 1 


yo+V 

0 

0 


“o 

0 

0 


- s o 

- _ 


fo 



L u 2 


( 3 -) 



and the scalar error equation is 

e l +k 2 e l +k l®l +k 0 e i = ^ a 0 + “l + °2^ x l + ( a 1 4 ' 2< V*l 


+ - Bjf - ( B 0 + 8 1 )r + u 1 + u c 


(4*) 


where 


k o “ a 0 + g 0 
k r a i + g i 
k 2 “ a 2 + *2 
3 2 

Let s +k 2 s +k^s+kjj = (s+X^) (s+X^) (s+Xj). 


Then 


(p + X^Xp + (p + X^e^ = 

(P+X 2 ) Cp+X 3 ) (<> 0 VVl +l<l 2 V 2 _ V3 _ V4 J 


Vo ' V “ l > 2V ( X 2 + X 3 ) * 2 V 2 


" dF ( *2 v 2 ) ~Vl - ( W Vl 


Where are unknown constants depending upon 
initial conditions. Dividing by (s+X«)(s+X-> 
yields 1 3 

2 4 

Wi“ £ Vi ~ £ 

1=0 i=3 


+ <(' 1 exp[-X 2 t] + << 2 expI-Xjt) 


( 8 *) 


Consequently, 


«. ** - (— ) v e 

i d 1+1 i 1 


(— ) v.e 

i+1 1 1 


The observer has the form 


1 - 0 , 1,2 


i - 3,4 


( 10 *) 


- k 2 

i 

o 1 


8 2 

0 

o’ 

- k l 

0 

1 

Z + 

8 1 

0 

0 

; k o 

0 

0 


. g o 

0 

0 


d^W + ( W *3 V 3 + *3*3 


+ dI ( ^3 V 3 ) + (X 2 + V *4 V 4 + *4*4 


+ d?< ^4 V 4^ + *1 + u 0 


'o ' 


0 ' 

b. 

r + 

u. 

1 


1 

b 0 


> 


(5*) 


where 


*0 “ a 2 


a 0 - X 2 X 3°2 


v. « X 
0 1 

V l + ^ X 2 +X 3^1 +X 2 X 3 V l"*l 


*2 “°1 " (X 2 +X 3 )a 2 y (X 2 +X 3^2 +X 2 X 3 V 2 “ *1 

<f 3 “Bj V 3 + ^ X 2 +X 3^3 +X 2 X 3 V 3 = r 

v 4 + < X 2 + X 3 ) V X 2 A 3 V 4‘ sr 


\ =B o 


Let 


»4 V 4 - *3 V 4 + Vl + *2 V 2 


u 0 ° V “ V 0 + * 2 + * X 2 + W + ^ i ^ i + ^ 2 + X 3 ^ v 1 ^ 

- « 3 fr + v 3 +( X 2 + X 3 ) v 3 ] - V V * A 2 + W 


Then 


(s+X^) (s+X 2 > (s+X 3 )e 3 = 


(s+A )(s+X ) [ l $ v - l tv 1 
*■ i *■ 1=0 1=3 J 

+ f n / 


(7*) 


where 

«2 

«1 

*0 

b . 


V" 


and 


n> 3 e l V 0 


- v 2 + (X 2 +X 3> V 0 )e l 

" V 1 + X 2 X 3 V e l 


— e- v - 
m 4 1 3 


— e i v / 

1 A 


1 0 
t 

a 2 -[ d 2 dt 1 


t c 

*1- 1 “,dt a - I d 

1 A * £ ft 


d 2 dt 1 


z (*) 



A Simulation 


The third order system of the example was 
simulated on a digital computer using the follow- 
ing parameters 


V 24 

P 

O 

H 

O 

Cj=30 

V 24 

mg/m^SOOO 

a 1 c 26 

a 1 =74 

C 2 =195 

k 1 =26 

mg/m ^.=200 0 

a 2 =9 

a 2 =0 

b 1 =30 

k 2 =9 

g 0 -g 2 =° 
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The eigenvalues of the observer (determined by 
{k^} were X^= -4, X 2 = -2 » * 3 “ The input to 

the plant was a square wave of magnitude 1 and 
frequency 6 t. Two parameters, b^ and g^ were 

adjusted by the adaptive law. These were initially 
set at bp = 73, g^ *=* -5 corresponding to a correct 
value of bg - 75, g^>= -74. The accompanying graph 

Illustrates the behavior of b^, g^, e 2> and e^ 

as a function of time. 


Remark 


A 'V 

As has been previously stated, m “ Tz and |im 
u - o). In the general case of an arbitrary plant 
matrix 3 l, the determinant of ^ may vanish for some 
instances of time. These momentary occurrences, of 
course, have no detrimental effect on m since con- 
vergence of u to id is guaranteed. In the important 
particular case of the preceding example, however, 
advantage has been taken of the fact that det T is 
constant by writing equation (*) as £ „ z. 

Since for the case of phase variable plant of high 


order the literal form of T 1 is easily produced, 

it is surmised that writing 'S ^ = T allows a 
particularly simple construction of u> when digital 
computation, rather than analog, is desired. 
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Reference Adaptive Control Systems with Feed- 
back and Prefilter Adjustable Gains," 4th 
Princeton Congress, March, 1970. 
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Proc. IEEE, Vol. 114, pp. 299-302, Feb. 1967. 

8 . Gilbart, J. W., and Monopoli, R. V., "A 

Modified Liapunov Design for Model Reference 
Adaptive Control Systems," Seventh Annual 
Allerton Conference, October, 1969. 
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and P. Kudva, "Adaptive Control Using Lyapunov* 
Direct Method", Yale Univ. Tech. Report No. 
CT-43, New Haven, Conn., "ct. 1971. 


10. Lion, Paul Michel, "Rapid Identification of 
Linear and Nonlinear Systems," AIAA Journal, 
Vol. 5, pp. 1835-1842, Oct. 1967. 


11. Carroll, R. L. "A Theorem About Constrained 
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Tech. Rept . 72-5, Dept, of Elect. Engineering, 
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Conclusion 

An adaptive observer has been demonstrated 
for single-input single-output systems with con- 
stant or slowly varying parameters. Work is 
currently underway to extend the observer to 
multivariable systems as well as Bystems with 
rapidly varying parameters. It is hoped that the 
adaptive ovserver will be eventually used not 
only for observing the state of an unknown system 
but in model reference problems and pole place- 
ment problems as well. 
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